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Abstract. Recently, it has been shown that it is possible to re- 
construct the projected mass distribution of a cluster from weak 
lensing provided that both the geometry of the universe and the 
probability distribution of galaxy redshifts are known; actually, 
when additional photometric data are taken to be available, the 
galaxy redshift distribution could be determined jointly with 
the cluster mass from the weak lensing analysis. In this paper 
we develop, in the spirit of a "thought experiment," a method 
to constrain the geometry of the universe from weak lensing, 
provided that the redshifts of the source galaxies are measured. 
The quantitative limits and merits of the method are discussed 
analytically and with a set of simulations, in relation to point 
estimation, interval estimation, and test of hypotheses for ho- 
mogeneous Friedmann-Lemaitre models. The constraints turn 
out to be significant when a few thousand source galaxies are 
used. 

Key words: cosmology: observations - gravitational lensing 
- dark matter - galaxies: clustering - galaxies: distances and 
redshifts 



1. Introduction 

It has long been recognized (see Tyson et al. 1984) that the 
gravitational field of a cluster acts on the images of distant 
galaxies by changing their orientation so that their major axes 
tend to become perpendicular to the direction of the center of 
the cluster. In turn (see Kaiser & Squires 1993) the mean ellip- 
ticity of the galaxy images can be used to measure the shear of 
the lens, a quantity related to the two-dimensional mass distri- 
bution (projected along the line of sight). By now a number of 
reconstruction techniques have been proposed and tested (see, 
e.g.. Kaiser et al. 1995; Seitz & Schneider 1996; Hoekstra et 
al. 1998). 

In the reconstruction process, information on the ratio of 
the distances from the observer to the source galaxies and 
from the lens to the galaxies is used. In practice, in the limit 
where the galaxies are much farther than the lens, the recon- 
struction process becomes independent of the distance of the 
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source galaxies (Schneider & Seitz 1995, Seitz & Schneider 
1995). This, of course, is very useful, as usually the redshifts 
of the sources are not available. If this limit cannot be taken 
(as it happens if the cluster itself has a relatively high redshift, 
e.g., Zd > 0.5), then one has to take into account the different 
redshifts of the background galaxies. The result of the gravita- 
tional lensing does depend on the geometry of the universe. 

Seitz & Schneider (1997) have shown that it is possible to 
reconstruct the mass distribution of the lens provided that both 
the probability distribution of galaxy redshifts and the geom- 
etry of the universe are known (e.g., the values of Vl and 17a 
in the standard Friedmann-Lemaitre cosmology). If other data 
(in particular, the angular sizes and luminosities of the source 
galaxies) are taken to be available, then the problem of the joint 
determination of the galaxy redshift distribution and of the lens 
mass density, for an assumed geometry of the universe, can be 
solved using the so called "lens parallax method" (Bartelmann 
& Narayan 1995). This latter method is based on the decrease 
of the surface brightness of galaxies and on the corresponding 
increase of the lens strength with redshift. A third possibility, 
of determining the properties of the lens and the geometrical 
characteristics of the universe from an assumed knowledge of 
the distances to the source galaxies, is not natural because of 
the prohibitive demands posed by the measurement of a large 
number of redshifts. However, this last option is actually most 
natural in principle, given the fact that redshifts are the quanti- 
ties that can be best measured directly. 

In recent years the technique to obtain photometric red- 
shifts (Baum 1962) has been greatly improved and tested. A 
large number of photometric redshifts of very faint sources 
{rriB > 25) can be obtained without additional telescope time 
(with respect to imaging). For the Hubble Deep Field more 
than one thousand photometric redshifts have been estimated 
(Lanzetta et al. 1996). Moreover, the technique can be rather 
accurate (provided that a good set of filters is used). For ex- 
ample, in a first test on a sample of 27 galaxies in the Hub- 
ble Deep Field, by comparison with the spectroscopic redshifts 
now available, more than 68% of the photometric redshifts have 
errors |Az| < 0.1, and all redshifts have errors |Az| < 0.3 
(Hogg et al. 1998). This suggests that in a not far future we 
might use photometric redshift information for galaxies lensed 
by a cluster. 
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In this paper we show how we can determine both the lens 
mass distribution and the geometry of the universe if the red- 
shifts of the source galaxies are known. The method outlined 
here could be used, in a Friedmann-Lemaitre universe, to get 
information on the density parameter ft and on the cosmolog- 
ical constant parameter ft^, or could be used to test different 
cosmologies. At this stage the study proposed here may be seen 
as a "thought experiment" because of the very high number of 
redshifts that are found to be required in order to constrain the 
geometry significantly. We will show below that the galaxies 
that one should observe can hardly be found behind a single 
cluster at the magnitudes currently attainable. However, it is 
possible to combine data from different clusters. Combining 
the study of 5-10 clusters would lead to statistically mean- 
ingful constraints. With such a device, the application of the 
present method might even become feasible in a not too far 
future. Surprisingly, simulations show that the method applied 
to observations considered to be realistic for current plans of 
the Next Generation Space Telescope (NGST) should provide 
significant constraints on the cosmological parameters even if 
based on a single cluster 

The paper is organized as follows. In Sect. 2 we recall the 
main lensing equations generalized to sources at different dis- 
tances and introduce the "cosmological weight," a quantity that 
gives the strength of the lens scaled to the relevant distances in- 
volved. The reconstruction of the cluster mass distribution and 
of the cosmological weight is described in Sect. 3. Section 4 
addresses the issue of the invariance properties associated with 
the reconstruction analysis. Section 5 introduces the determina- 
tion of the cosmological parameters by recognizing the limits 
on the measurement of the cosmological weight. The results of 
a wide set of simulations are described in detail in Sect. 6. Sec- 
tion 7 addresses the problem of the actual feasibility of a mea- 
surement based on the procedure outlined in the paper Finally, 
in Sect. 8, we summarize the main results obtained. Four ap- 
pendices contain detailed discussions and derivations of some 
important results. 

2. Distance dependent lensing relations 

This section is aimed at defining the basic mathematical 
framework and mostly follows the article by Seitz & Schnei- 
der (1997). Let us consider a lens at redshift with two- 
dimensional projected mass distribution where is a 
two-dimensional vector representing a direction on the sky (the 
region of interest on the sky is very small and can be considered 
flat). For a source at redshift z we define the critical density 
^c{z) as 



oo for z < Zd 

Sc(^) — — — ^-J- r Otherwise . 

4:TTGD{zd)D{zd,z) 



(1) 



Here D{z) and D{zd,z) are, respectively, the angular 
diameter-distance from the observer to an object at redshift z, 
and from the lens to the same object. The quantities D{za), 
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Fig. 1. The weight function for a lens at Zd = 0.8 in four dif- 
ferent universes. From the top to the bottom: = 0.3, JIa = 
(dashed); ^ = 0.3, Vl^ = 0.7 (dotted); = 1, f^A = (solid); 
= 1, r^A = 1 (long dashed). 



D(z\ and Diz^Tz) replace the Euclidean distances, respec- 
tively Dqa, Dos, and D^s, used in Paper I (Lombardi & Bertin 
1998a). 

It is useful to define a redshift independent critical density 

as 



= lim Ee(^) , 



(2) 



related to the redshift dependent critical density through a "cos- 
mological weight" function 



w{z) — 



I]c(z) D{z) 



(3) 



In principle, the above relations could be the starting point also 
for investigations based on non-standard cosmological models, 
using the constraints on the cosmological weight w{z) pro- 
vided by the gravitational lensing analysis to be described be- 
low. In practice, for the rest of the paper we will focus on the 
standard Friedmann-Lemaitre cosmology, where the angular 
diameter-distance D{z) can be written as (Peebles 1993; see 
also Kayser et al. 1997) 



D{z) = D{0,z) 
D{zd,z) 



Ho{l + z)Vih 



: sinh 




where 



E{z') ^ ^n{i + z'f + nR{i + z')^ + nA , 



(4) 

(5) 

(6) 



and, we recall, il+ilA+rifl — 1. The hyperbolic sine in Eq. (||) 
is characteristic of an open universe (flu > 0). 

In general, w{z) — > 1 for z ^ oo. As we will see shortly, 
the lensing equations depend on the geometry via the cosmo- 
logical weight, and thus this is the function that can be deter- 
mined from lensing observations. The function depends in a 
complicated manner on the two cosmological parameters O and 
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f2A- For an Einstein-de Sitter universe (fl — 1, f^A — 0, so that 
— 0) one has, for z > z^, 



HoD{zd,z) 



1 



1 + 2 V\/TT^ VTTz 



w{z) 



x/m- 1 



(7) 



(8) 



Some examples of weight functions for a lens at redshift z^ — 
0.8 are given in Fig. I. It should be noted that the cosmological 
weights for "reasonable" universes lie very close to each other 
This of course makes the determination of the cosmological 
parameters from the measurement of the cosmological weight 
a difficult task. 

For sources at redshift z, we define the Jacobian matrix of 
the ray-tracing function as (see Eq. (3) of Paper I) 



1 



- 7i) 

W{z)'-f2 



W{z)j2 

1 - ■w{z){k + 7i) 



(9) 



Here k = k(0) = E(0)/S]^ is the lens density in units of the 
redshift independent critical density defined in Eq. (^), while 
7 = jid) is the complex shear. Then the redshift dependent 
amplification ij,{9,z) is given by 

|dctA(0,z)| 

li[6,z) 



[i-w{z)K{e)] -\ w{z)j{6)\ 



(10) 



This expression gives the magnification produced by the lens 
on a source located at redshift z. If this quantity can be mea- 
sured (for example by using galaxy counts), the mass density of 
the lens can be constrained strongly, and an otherwise present 
invariance property (see Sect. 4) can be broken. 

For a galaxy at redshift z seen at position 9 a complex quan- 
tity e is defined in terms of the galaxy quadrupole moment Qij 
(see Schneider & Seitz 1995 for a detailed discussion): 



Qii - Q22 + 2iQi2 



Qu + Q22 + 2v/QiiQ22-Q 



(11) 



At variance with Paper I and Paper II (Lombardi & Bertin 
1998b), it is convenient to measure the ellipticity using this 
quantity instead of x (see Eq. (7) of Paper I), because of the 
simpler transformation properties of e (see Eq. (|l4|) below). In 
fact, as shown by Seitz & Schneider (1997), the observed ellip- 
ticity is related to the source ellipticity through the relation (the 
superscript stands for "source") 



l-g*{e,z)e' 

i~g{e,zy* 



g*ie,z) 



for \gie,z) \ < 1 



otherwise 



(12) 



where the reduced shear 5(0, z) is defined as 

w{zhie) 



9(0, z) 



1 - w{z)K{e) 



(13) 



From the observation of a large number of galaxies near 9, 
with redshifts close to z, we can thus determine the reduced 
shear g{9, z). In fact, it can be shown that the expected mean 
value of the observed ellipticities is given by 



{e){9,z) = 




if |g(0,z)| < 1 
otherwise , 



(14) 



under the isotropy hypothesis (e*^) = for the source galaxy 
distribution. This expression depends on the redshift of the 
galaxies used. Therefore, if we measure the mean value of the 
ellipticities for galaxies at different redshifts, we can obtain in- 
formation on the cosmological weight function w(z). Once the 
cosmological weight function w{z) has been determined, the 
cluster mass reconstruction can be performed using an iterative 
procedure on Eq. (|l3|), seeded by k = 0, and some kernel anal- 
ysis (see Paper II). As shown by Seitz & Schneider (1997), this 
can be carried out even if one does not have the redshift of each 
galaxy, provided one knows the galaxy probability distribution 
in z. 

In the rest of the paper we assume that the ellipticity disper- 
sion of the source galaxies is independent of redshift (see also 
App. C). 



3. Joint determination of the shear and of the cosmological 
weight 

Let us now address the goal of this paper, the joint determi- 
nation of the lens mass and of the geometry of the universe. 
In the following, for simplicity, we will consider a cluster sub- 
critical for all z, i.e. a cluster that cannot produce multiple im- 
ages, characterized by \g{9, z)| < 1 for all 9 and all z (from 
Eq. (10), by recalling that subcritical lenses are also character- 
ized by det A > {), it can be shown that this happens when 
n{9) + |7(0)| < 1 for every 9). General clusters could be 
treated basically in the same way, but with considerably heavier 
notation. Let the positions 0*^"^ the observed ellipticities e^"^, 
and the redshifts z^"^ of N source galaxies (71 = 1, ... , TV) be 
known, together with the redshift Zd of the cluster 

In the following we will use subscript for true quantities. 
For example, ko will be the true dimensionless surface density 
of the lens, and Wo(z) the true cosmological weight (given by 
Eq. (|^)). For the complex quantities 7 and g we will use the 
notation 7oi and goi to denote the i-th component of the true 
value. Following the notation used in statistics, we will use hats 
for measured quantities ("estimators"). Thus k{9) will be the 
measured mass density, and i«(z) the measured cosmological 
weight. 

In order to describe the reconstruction process, we start by 
considering the weak lensing limit, and then we generalize the 
results obtained. 
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3.1. Weak lensing limit 



3.2. General case 



The weak lensing limit is characterized by |7o(0)|, ko(0) ^ 1, 
so that the reduced shear .go(^i z) can be identified with the 
product wq{z)^q{6). In this limit Eqs. ( [l2[ ) and ( p4| ) become 



lo[0)wo{z) 



and 



(e)(0,z) = -7o(0)u;o(2:) 



(15) 



(16) 



(17) 



Here, we recall, the cosmological weight function is taken to 
be known. 

Similarly, if we take the shear 7(0) to be available, we can 
introduce a redshift weight function Wz{z, z'), with the condi- 
tion that Wz{z, z') is significantly different from zero only for 
\z — z'\ small. Then the cosmological weight function w{z) can 
be estimated by 

i:n=iW.{z,z(^))^[^{e^^'>)e(^)*] 



'w{z) = — - 



(18) 



where and the asterisk denote real part and complex conju- 
gate (we recall that e and 7 are complex quantities). 

The two equations can be used to obtain by iteration the 
joint determination of the shear and of the cosmological weight 
function. We start with a guess for the weight function w{z) 
(e.g., the expression known for the Einstein-de Sitter cosmol- 
ogy given in Eq. (^), and calculate the shear using Eq. (|l7|). 
Then we insert the estimated shear 7(0) in Eq. (|l8|) and thus 
obtain a first estimate for 'w{z). We then go back to Eq. ( p^ 
and obtain a new determination of 7(6), etc. 

As a by-product of the analysis, by means of standard meth- 
ods (Seitz & Schneider 1996; Paper II), we can also recon- 
struct, from the estimate 7(0), the dimensionless mass distri- 
bution 

At the end of the iteration process, it can be shown that 
the procedure leads to the correct determination of w{z), in the 
sense that 



{w){z) 



z, z 



{ny 



Wo (2 



i:L,Wz{z,z(^^)) 



Wq{z) . 



(19) 



Further discussion of Eqs. ( p^ and (|l8|) is given in App. C. 



If the lens is not weak the solution of our problem is more com- 
plicated and requires nesting two different iteration processes 
Let us start by assuming that the cosmological weight is known 
A relation between 7(0) and the observed ellipticities now in 
volves the quantity k{6), as can be easily found from Eqs. 
andQ: 



N 



These relations allow us to estimate, from a set of measure- 
ments one of the two quantities, either the shear or 7(6) 
the cosmological weight function, provided the other is known. 
The basic idea is to average the observed ellipticities of galax- 
ies close to each other either in 0-space or in redshift, following 
the general theory of maximum likelihood methods to mini- 
mize errors (e.g., see Eadie et al., 1971). 

If we take w{z) to be known, we can introduce the spa- 
tial weight function W{6, 6') with the property that W{0, 6') 
is significantly different from zero only for \\6 — 6'\\ small 
(see Kaiser & Squires 1993; Paper II). Then, from Eq. ([l5|), 
the shear 7(0) can be estimated by 



w 



(z(")) 



1 - 



,(") 



N 



n=l 



1 - 



(20) 



It is not difficult to recognize that this equation, which gener- 
alizes Eq. ([TtI), leads to the optimal estimation of 79. We note 
that the combination of weights is again optimized with regard 
to the error on 7. It is clear that Eq. ( pO| ) can be used to ob- 
tain both the shear 7 and the mass density n, provided that we 
know the cosmological weight. To this purpose we can use an 
initial guess for k{6) (e.g., = 0) in Eq. (^0|), obtaining 
in this way a first estimation of 7(6). From this we can ob- 
tain, using standard methods (see Kaiser 1995, Seitz & Schnei- 
der 1996, and Paper I and II for a discussion), the correspond- 
ing mass density, up to a constant. Then we can go back to 
Eq. ( pO| ) with the mass density and obtain a new determination 
of 7, and so on. The method outlined here is similar to that 
described by Seitz & Schneider (1997), with the important dif- 
ference that they cannot use the weights in Eq. ( pO| ) 
since the galaxy redshifts are not taken to be known. 

For given values of 7(0) and k{0) we can estimate w{z) in 
a manner similar to the weak lensing limit. In the general case 
Eq. dl^ ) is replaced by 



w{z) 



N 



Wz{z 



[1 



w[z 



r N 



7(0(")) 



r,(0(")) 



(21) 



As in the weak lensing limit, the joint determination of 7 
and w is performed by an iterative process between Eq. ( pO[ ) 
and Eq. (pi)). However, we stress once more that in this case 
Eq. ( pc| ) requires a separate iteration for the joint determina- 
tion of 7 and n, to be performed between Eq. ( pO| ) and some 
kernel analysis for mass reconstruction (see Paper II). 

The proof that Eq. ( pi] ) leads to an estimation of 'Wq{z) 
would be similar to that for Eq. ( 18 ). 



4. Invariance properties 

So far, we have assumed that the shear map 7(0), the mass 
distribution k{0), and the cosmological weight w{z) can be 
determined uniquely. Unfortunately, two invariance properties 
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affect the lensing equations and leave the reconstruction pro- 
cess ambiguous. 

The first property is the well known "sheet-invariance" 
(Kaiser & Squires 1993; Schneider & Seitz 1995). In the weak 
lensing limit, this affects only the dimensionless mass density 
k{9) and has a simple interpretation. If we add a homogeneous 
"sheet" to the mass density all equations remain unchanged. 

The second invariance property, that we call "global scal- 
ing invariance," is a new feature associated with the process 
of the joint determination of the lens mass density and of the 
cosmological weight. 

4.1. Weak lensing limit 

In this limit the two invariance properties refer to the transfor- 
mations (sheet) 

k(6/) ^ K{e) + C , 7(6») ^ 7(6») , w{z) ^ , (22) 
and (global scaling) 

k(9) kK(e) , 7(0) ^ k-f(e) , w(z) ^ . (23) 

k 

We note that, in the weak lensing limit, the two invariance prop- 
erties are decoupled, in the sense that the first involves only the 
dimensionless density and the second is met in the determina- 
tion of the cosmological weight (see Sect. 3.1) even without 
considering the estimation of the dimensionless density. 

4.2. General case 

In the general case, the two invariance properties are coupled, 
because all three basic quantities k, 7, and w have to be deter- 
mined jointly. In particular, when the lens is not weak, a new 
undesired source of uncertainty in the determination of the cos- 
mological weight is introduced by the sheet invariance itself. 

The sheet invariance arises because of the differential re- 
lation between 7(6) and k{9). Let us recall the method used 
to obtain the dimensionless density map k{6) from the shear 
7(0). First the vector field (Kaiser 1995) 

u{9)^-(^'''+^^A (24) 

V">'2,1 - 71,2 / 

is calculated from the shear map. As shown by a simple calcu- 
lation, this vector field is the gradient of the mass density, i.e. 
'Vk{6) = u{6). Hence, from the shear 7, the density k can be 
determined only up to a constant. This is the usual sheet invari- 
ance in the weak lensing limit, where the shear is the observed 
quantity. 

In the general case one cannot measure 7(0) directly be- 
cause Eq. (pO|), used to obtain the shear, involves the dimen- 
sionless mass density. In addition, a simple generalization to a 
relation Vk = u (see Kaiser 1995 and Paper II) is not avail- 
able because of the presence of the cosmological weight (see, 
however, Seitz & Schneider 1997 for an approximate relation). 
Thus an iterative process is needed. In the iteration, the arbi- 
trary constant used to invert the vector field u into the density 



map K and the cosmological weight act in a complicated man- 
ner, which is hard to describe analytically. If the lens is sub- 
critical and if the cosmological weight is left unchanged (i.e. 
if the outer iteration is not performed, see Sect. 3), the sheet 
invariance can be expressed by (see Appendix A) 

K^^j^ + {l-C)Kie), 7(0) ^ (1 - C)7(0) , (25) 

where {■)z is the average over the redshifts. A similar expres- 
sion has been found by Seitz & Schneider (1997), by consid- 
ering the sheet invariance in a different context (the lens re- 
construction for a population of galaxies with known redshift 
distribution). Unfortunately, even under the above special con- 
ditions, if one introduces the transformation ( p5| ) in Eq. (pT]), 
the resulting measured cosmological weight w{z) changes in 
a complicated manner. It can be shown that the related trans- 
formation for the cosmological weight cannot be reduced to a 
simple scaling w{z) t-^ vw{z). 

The global scaling invariance property is a the same as that 
encountered in the weak lensing limit; but here we recall that 
the three quantities 7, w, and n have to be determined jointly. 
In particular, if we consider the transformation 

k(6») ^ kK(e) , 7(0) ^ k-f(e) , w(z) ^ , (26) 

k 

all equations remain unchanged. We observe that a similar in- 
variance property is met when one derives information on the 
geometry of the universe from the separation of multiple im- 
ages produced by a strong lens. 

4.3. General strategy 

In conclusion, there are two invariance properties, sheet and 
global scaling invariance. In the general case, both invariance 
properties affect the determination of the cosmological weight 
■w{z): the global scaling acts as a simple scaling, while the sheet 
invariance acts in a different and a priori unknown manner. 
Therefore, to the extent that the impact of the sheet invariance 
remains out of control, any determination of the cosmological 
weight would be useless. This implies that when the lens is not 
weak we can choose between two possibilities. Either we dis- 
card the central part of the cluster, leaving only the parts where 
the lens is weak, or we must find a way to break the sheet in- 
variance. 

As discussed by Seitz & Schneider (1997), there are several 
possible ways to break the sheet invariance. A simple and use- 
ful method is based on the magnification effect (Broadhurst et 
al. 1995), i.e. on the observed (local or total) number of galax- 
ies. In the case of a single redshift source population, a simple 
expression holds (see Appendix B) for the total projected mass 
associated with k 

M = J k{0) d^e* = i y [1 - det A{e)\ d^e . (27) 

This has a very simple interpretation. Suppose that galaxies are 
distributed on the source plane with constant number density 
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p^. Then, the observed number density p{6) (number of galax- 
ies per square arcsec) is given by 



p{e) = p^ldct A(6>)| = 



p{e) 



(28) 



Inserted in Eq. ( |27| ) it gives 

M = ^ 1 [p'-p{e)sguA{e)\d' 



(29) 



In this expression sgn A{9) — sgn(det ^(0)) is the sign of 
the determinant of A{6), Ni is the observed number of images 
of kind i (see Schneider et al. 1992 for the relevant classifica- 
tion of images), and iVcxp the expected number of galaxies in 
an area equal to that used for the observations. (We note that 
the classification of an observed image is easy once the critical 
lines can be determined empirically.) 

In view of the above discussion, in the following we will 
always consider either the weak lensing limit or the general 
case with the sheet mass degeneracy broken. Simulations show 
that, for the present problem, a (not too "strange") lens with 
K < 0.3 everywhere can be considered to be weak, i.e. to meet 
the asymptotic requirements k, I7I <C 1. The global scaling 
invariance, instead, will be addressed explicitly in the analysis, 
i.e. it will be considered present and not broken. 

One might think that the global scaling invariance could be 
broken simply by the use of the asymptotic behaviour of w{z) 
ai z ^ 00, i.e. by setting the cosmological weight equal to 
unity at high redshifts. Unfortunately, the asymptotic regime is 
reached only slowly. Thus only galaxies at very high redshifts 
(say z > 10) would be useful, which is practically beyond the 
limits of observations. Still, in the discussion of the simulations 
(see Sect. 6), we will show that the fact that the global scaling 
invariance, in contrast with the sheet invariance, can be reduced 
to a simple scaling, allows us to proceed to a meaningful deter- 
mination of the cosmological parameters (once the sheet invari- 
ance is assumed to be broken). In fact, we will fit the shape of 
the measured cosmological weight with a set of cosmological 
weights (corresponding to different cosmological parameters), 
leaving unspecified the scale of the function. 

5. Limits on the determination of the cosmological weight 
function 

In principle, given a determination of the cosmological weight 
function w{z) from the procedure outlined in Sect. 3, it is 
straightforward to derive information on the parameters and 
f^A that characterize the assumed Friedmann-Lemaitre cosmol- 
ogy. In particular, we may resort to a simple least squares or 
maximum likelihood method. Such a method is able to lead not 
only to "point estimations" of the cosmological parameters, but 
also to "interval estimations," based on the identification of the 
relevant "confidence level" regions (see Appendix D). In prac- 
tice, difficulties arise because the determination of w{z) will 
remain incomplete: 



- The scale of the cosmological weight is expected to be un- 
known, because of the global scaling invariance described 
in Sect. 4. 

- The measurement of ^(2) is going to be restricted to a lim- 
ited redshift interval and we anticipate to be able to do so 
with only limited redshift resolution. 

- The measurement of w{z), resulting from the study of a 
finite number of source galaxies, is going to be affected by 
statistical errors. 

Therefore, we now need to determine, for every cosmologi- 
cal model, both the expected measured cosmological weight 
(w)(z) and the related error (covariance matrix). As we de- 
scribed in Paper II, since we are dealing with functions, ei- 
ther of space k{0) or of redshift w{z), the covariance ma- 
trix analysis has to be extended to the related two-point cor- 
relation functions. These calculations can be carried out either 
analytically (see Appendix C) or numerically, through Monte 
Carlo simulations. In any case, we may anticipate that while 
an error on the expected measured {w){z) introduces a bias 
in the estimation of the cosmological parameters, an error on 
the (co)variance of the cosmological weight is only bound to 
increase the expected variance of the cosmological parameters 
thus determined, which is not critical for point estimations. 

6. Simulations 

6.1. Method 

We refer to a cluster at redshift Zd = 0.8. Different geome- 
tries of the universe set different values for the critical den- 
sity S^. In particular, calling h = i?o/(100 km s"^ Mpc"^) 
the reduced Hubble constant, the case = 1, Ha — has 
Y.^ = 5.49 h kg m^^ while the case n = 0.3, f^A = 0.7 has 
— 3.99 h kg m~^. In the following we specify the cluster 
mass distribution directly in dimensionless form. For this we 
consider the sum of a diffuse component A (a truncated Hubble 
profile) and three compact objects B, C, D. The corresponding 
(projected) two-dimensional distributions are 



- ) .2 2 for r < i? , 

otherwise , 



and 



Ko 



1/ro - 1/R \ ^rl + r2 ^fW~W^ 



(30) 



(31) 



Here r denotes the angular distance from the center of each 
object. For simplicity of notation, we use the same symbols 
Kq, ro, R for three scales, taken to be different in each case 
(see Table 1). The dimensionless density n(S) is calculated on 
a grid of 50 x 50 points (corresponding to a field of 5' x 5' on 
the sky). The cluster, shown in Fig. 2, has a mass (for h = 0.5, 
= 1, and r^A = 0) within a 1 Mpc aperture diameter of about 
2 X lO'^^Mo. From the lens equations we calculate the shear 
7 on the same grid. Within the Friedmann-Lemaitre geometry. 
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the weight function is calculated directly from Eqs. (3), (4), and 
(5). For the purpose, 25 equally spaced points in redshift from 
Zd = 0.8 to Zd = 5.8 are found to be sufficient. 

The ellipticities of source galaxies are generated using a 
truncated Gaussian distribution 



1 



7rP2[i_exp(-l/P2 



■ exp 



,s|2 



(32) 



with P = 0.15. Note that the distribution is normalized with 
the condition 27r /p(|e''|) le**! dje*^! = 1. Galaxy positions are 
chosen randomly on the field. For simplicity, following Seitz 
& Schneider (1996), we take a uniform distribution in the ob- 
served position, thus neglecting magnification effects. Finally, 
following Brainerd et al. (1996), we draw the redshifts from a 
gamma distribution 



Pz{z) 



2zl 



exp(-z/zo) 



(33) 



where zo = 2/3. This distribution has maximum at z = 2zo ~ 
4/3 and mean (z) = Szq = 2. Current surveys favor a slightly 
higher value of zq (Driver et al. 1998, for a limiting magnitude 
/ < 26). All distributions can be obtained from a uniform dis- 
tribution by the transformation or rejection methods (see Press 
etal. 1992). 

Then the observed galaxy ellipticities are calculated using 
Eq. (p^. The "observations" used for the reconstruction pro- 
cess are then the set of galaxy positions {0^"^ }, the set of ob- 
served ellipticities {e'") }, and the set of redshifts {z^")}. The 
simulations are carried out with Gaussian weights in ||0 — 0'| 
(with scale lenght 9") and in |z — z'| (with scale 0.2). 

The reconstruction process follows the method described in 
Sect. 3, and can be summarized in 9 steps: 

1. To begin the reconstruction, an Einstein-de Sitter uni- 
verse is assumed, and thus, as first guess, the cosmological 
weight (z) is calculated (Eq. ^). 

2. The initial density kI^'"! and the initial shear ^[O'"] are both 
set to zero. 

3. A first determination of the shear --ylO'il {6) is then obtained 
from the observed ellipticities using Eq. (pO[). 

4. The new mass distribution k}-^'^\6) is calculated using the 
shear 7[O'il(0). The inversion from the shear into the mass 
distribution is performed by convolution with the curl-free 

Table 1. The parameters used for the four components of the 
cluster mass distribution. Component A is a Hubble profile, 
components B, C, and D represent compact objects ("domi- 
nant" galaxies). 



kernel H^^{e, 6') (Seitz & Schneider 1996; for a discus- 
sion on the merits of different kernels see Paper II). The 
sheet invariance is broken by requiring the total recon- 
structed mass J k{0) d?6 to be equal to the true value (see 
Sect. 4.3). 

5. Points (2), (3), and (4) are iterated. This is repeated a num- 
ber / of times sufficient to obtain good convergence (typi- 
cally I — b for sub-critical clusters). 

6. The final determinations kI"'^! and '^[O-^l are used in 
Eq. ( pi] ) to get the new weight function w^^l . 

7. Points (2-5) are iterated, thus leading to the determination 
of ixil^+^l from the knowledge of w'^l. (The outer iteration 
on i involves w, while the inner iteration on i involves n 
and 7.) As starting density and shear for the new inner iter- 
ation we use the last determinations of these quantities, i.e. 
^[j+1,0] ^ ^jjjj simiiai-iy ^[j+i.o] = ^[j./], xhe outer 

iteration is performed a number J of times. 

8. The final lens density ^[^■^1 is then compared to the true 
one kq. 

9. The final cosmological weight wl'^l is then used to deter- 
mine, or at least constrain, the geometry of the universe. 
The methods used are the ones described in detail in Ap- 
pendix D. 

Simulations show that, for "reasonable" values of J7 and 
r^A, the factor k related to the global scaling invariance remains 
very close to unity; thus the measured cosmological weight 
does not differ significantly from the true cosmological weight. 
This happens because the cosmological weight changes only 
slightly for different values of Q, and J^a (cfr. Figure 1). Hence 
the initial choice w[°l (z) (Einsten-de Sitter universe) is not far 
from the true wq (z) . In this situation it is correct to use the mea- 
sured total mass to break the sheet invariance, as done in point 
4 (instead, in the presence of a wrong scaling factor fc, this pro- 
cedure might lead to undesired results because of the coupling 
of the sheet invariance with the global scaling invariance). 

6.2. Results 

A wide set of simulations show that the method produces sig- 
nificant results if the number of galaxies studied is a few thou- 
sand or more. 

The first part of the reconstruction deals with the shear ^{0) 
and with the mass distribution k{9). Figure 2 shows an exam- 
ple of mass density obtained from the method. The map ob- 
tained is rather detailed and smooth because of the large num- 
ber of galaxies used {N = 20 000). As explained in detail in 
Paper II, these properties depend on the spatial weight function 
W{0, 6') used: the larger the characteristic size of this func- 
_tion, the smoother the shear and density maps obtained. This 
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Density behavior is consistent also with the calculations of Appendix C. 
Jta the simulations this size has been adjusted to the value that 
gives the best results on the following discussion of the cosmo- 
logical parameters. 

The initial guess for the cosmological weight function ui'^l 
js not very important for the determination of the shear, because 
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Fig. 3. The weight expected mean and error on the weight 
function for an Einstein-de Sitter universe. The solid line 
shows {w){z), the dashed lines the expected errors (obtained 
as •\/Cov(w; z, z)) and the dotted line the true wo{z). The er- 
rors are calculated for = 10 000 galaxies; errors for different 
values of N can be obtained by recalling that the error scales 
as %/iV. 



Fig. 2. The original dimensionless mass density and an exam- 
ple of reconstruction. The number of galaxies used is iV = 
20 000; the assumed cosmological model is Einstein-de Sitter 
The cluster occupies a square with side of 5', corresponding 
approximately to 1.23 h^^ Mpc. 

a wrong choice would only lead to an increase of the errors (at 
least in the weak lensing limit; see comment after Eq. (C.4)). 
This behavior is well verified in the simulations. At the end of 
the first inner iteration, the mass density obtained, k^^'^\ does 
not differ significantly from the determinations obtained at the 
end of the following cycles, kI-^'^I (except for a factor arising 
from the global scaling invariance). After the first inner cycle, 
the mass and the shear maps are close to the best ones that 
we can hope to have. This suggests that the number of inner 
iterations / could decrease after the first cycle. Thus we could 
start with 1 = 5 for j = 0, and then let / = 3 or even 1 = 2 for 
j > 0. The reduction of / can lead to a significant reduction of 
machine-time. 

The following step is the determination of the cosmolog- 
ical weight. Simulations show that the number of outer itera- 
tions needed to obtain a good estimation of the cosmological 
weight is very low, say J = 3, because the method is able to 
give the correct shear map after the first inner loop. Simulations 
also show that this property, strictly expected only in the weak 
lensing limit, is in fact valid (at least approximately) in the gen- 
eral case, provided the lens is not too "strong." The estimated 
cosmological weight w(z) (solid lines in Fig. 4) is smooth on 
the characteristic scale of Wz {z, z') (see Eq. (C.7)); because of 
this, it remains at a finite value at z = z^. Moreover, as ex- 
pected, the error on the cosmological weight increases slightly 
for z near z^ and more for high values of z, where the number 
of galaxies decreases (see Fig. 3). In particular, the smooth- 
ing of the cosmological weight is very important near the clus- 



ter, at z ~ 2d, where the true cosmological weight vanishes 
abruptly. This clearly indicates that neither the limit w{z) 1 



for z — > oo nor the limit for z 



can be used to break the 



global scaling invariance. Figure 4 shows the reconstruction of 
the cosmological weight in the case of an Einstein-de Sitter 
universe; different choices for il and ft^ lead to similar results. 
Figure 4 should be compared with Fig. 3 which shows the ex- 
pected mean value (w) (z) and the related error for an Einstein- 
de Sitter universe with = 10 000 galaxies. Figure 3 clarifies 
the general properties discussed above and, in more detail, in 
Appendix D, for the expected error on w{z) and suggests that 
our method should be able to constrain significantly the cosmo- 
logical parameters. In fact, the error on w{z) for z ~ 2 is suf- 
ficiently small to distinguish different consmological models 
even if the measured weight are affected by the global scaling 
invariance. 

From the estimation w{z) of the cosmological weight we 
can obtain information on the cosmological parameters as ex- 
plained in Appendix D. For a description of th e results, we plot 
the contours of the £2 function of Eq. ( ^.^ ) in a square do- 
main [0, 1] X [0, 1] of the r^-ilA plane. Figures 5 and 6 show 
the confidence regions obtained for various confidence levels 
CL in different cosmological models and with different num- 
bers of source galaxies. From diagrams of this type we argue 
that = 10 000 can be considered as a lower bound for the 
applicability of our method. 

These figures also show that, unless an exceedingly high 
number of source galaxies is available, point estimation is not 
very meaningful. In fact, these examples show that the mini- 
mum of the function can occur quite far from the true values 
of the cosmological parameters. On the other hand, for not un- 
reasonable (see discussion in Sect. 7) values of N, the method 
does constrain the cosmological parameters. If additional a pri- 
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Fig. 4. The reconstructed weight function in an Einstein-de Sit- 
ter universe (17 = 1, 17 a = 0). The dotted line represent the 
true weight function w^iz); the dashed line is the expected 
mean value as predicted by Eq. ( |C.7| ); the solid line is 

the measured The three figures have been obtained with 
= 10 000 (top), AT = 20 000 (middle) and TV = 50 000 
(bottom) of source galaxies. 



ori conditions are included (e.g., if the universe is taken to have 
= 0), parameters can be constrained rather efficiently. 
In order to check the distribution of we have repeated a 
large number of times the entire process, using each time a dif- 
ferent set of source galaxies. The measured probability distri- 
bution of I2 for different values of N is shown in Figs. 7 and 8. 
In these figures, the histograms show the number of simulations 
that have produced a value of £2 in the corresponding interval. 




Fig. 5. The confidence regions corresponding to CL — 68%, 
CL = 80%, CL = 90% and CL = 95%. The true value of the 
cosmological parameters, Vl — 0.3 and J^a = 0.7, is shown as 
a white cross. The galaxies used are = 10 000 (top), iV = 
20 000 (middle) and iV = 50 000 (bottom). 
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Fig. 6. As Fig. 5, but with Q = 1 and = 0. 



Fig. 7. The observed probabiUty distribution for ^2- The 15 bins 

have been chosen so that the expected number of events per bin 
would be constant if ^2 followed a chi-square distribution with 
one degree of freedom. The graphs show that, except for the 
first two bins, the observed probability distribution is roughly 
constant. The three graphs correspond to = 10 000 (top), 
TV = 20 000 (middle), and TV = 50 000 (bottom). 

Intervals are chosen so that the theoretical number of events 
per bin is constant. From these and other results not shown it 
is clear that £2 tends to a chi-square distribution with two de- 
grees of freedom only for a very large number of galaxies iV, 
say N > 100 000. In contrast, for "small" values of N the 
observed distribution is (approximately) consistent with a chi- 
square with one degree of freedom. This suggests that, when 
the number of galaxies is small, a new approximate invariance 
may be present. The method is able to break this invariance for 
larger values of N. 

This behavior is also suggested by the shapes of the confi- 
dence regions of Figs. 5 and 6. In fact, the confidence regions 
obtained have a characteristic shape approximately following 
the lines -|- = 1 — = const. Thus our method should 
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Fig. 8. As for Fig. 7, but with bins chosen so that the expected 
number of events per bin would be constant if £2 followed a 
chi-square distribution with two degrees of freedom. 

be able to constrain well the curvature flu, while the quantity 
rt — rt\ should remain practically unknown, at least when the 
number N of galaxies used is not too high. This approximate 
invariance may explain the behavior of the quantity £2 for finite 
N. In fact, if the invariance were exact, £2 would be expected 
to follow a chi-square distribution with (only) one degree of 
freedom. This is indeed observed for N < 50 000. Thus in or- 
der to draw the confidence regions, one should use the results 
of Monte Carlo simulations, unless the number N of galaxies 
is very large. The contours presented in this paper have been 
obtained from Monte Carlo simulations. The use of the asymp- 
totic limit would have led indeed to larger confidence regions. 

7. Feasibility 

In this section we briefly discuss the applicabiUty of our 

method to observations. As noted earlier, the main difficulty 
with our method is probably the high number of galaxies (and 



redshifts) needed to obtain significant constraints on the geom- 
etry of the universe. 

In order to show that the examples chosen for our simula- 
tions can be considered to be not unrealistic, we first note that 
the required density of galaxies is between 400 and 800 galax- 
ies per square arcminute. This density is certainly too high for 
current instruments. In fact, the current achievable density of 
galaxies (with "reasonable" integration times) is of the order 
of 60 galaxies per square arcminute (Clowe et al. 1998), and 
this number is only likely to double with the Advanced Cam- 
era on HST. However, current estimates for NGST (Stiavelli et 
al. 1997) give more than 1 000 galaxies per square arcminute. 
It is also interesting to note that the baseline for the camera cur- 
rently considered for NGST is about 4', not far from the value 
used in our simulations (5'). 

A second difficulty with the method proposed here is the 
large number of redshifts required. As suggested in the Intro- 
duction, photometric redshifts might be used. In this case, how- 
ever, sizable statistical errors on the measured redshifts are ex- 
pected. At this stage, for simplicity, such errors have not been 
taken into account in the simulations, but they obviously should 
be considered in detail before applications are tried. (We note 
also that errors on the measurement of ellipticities may have an 
impact, especially if small and weak source galaxies are used.) 

Another important ingredient is the cluster chosen in our 
simulations. In an Einstein-de Sitter universe with h = 0.5, 
the cluster would have a mass within a 1 Mpc aperture of about 
2 X lO'^^ Mq. We may compare our cluster with the X-ray clus- 
ter MS 1054 - 03 located at Zd = 0.83 (Luppino & Kaiser 
1997); weak lensing estimates for this cluster give a mass 
within 1 Mpc in the range 1.2 x 10^''' Mq and 5.5 x 10^^ M©, 
depending on the assumed redshift distribution for the source 
galaxies. 

In conclusion, the proposed method should give interest- 
ing results if applied to NGST observations. Of course, several 
refinements and improvements would be required for the pur- 
pose. For example, as suggested in the Introduction, the use of 
several clusters should provide tighter constraints on the cos- 
mological parameters. Broadly speaking, the use of A^ci similar 
clusters should make the errors on the cosmological parameters 
smaller by a factor \/Nci. However, for this purpose the method 
should be generalized to be applied to different clusters at dif- 
ferent redshifts (actually, we expect that if the redshifts of the 
clusters are different, more information on the geometry of the 
universe is available). We should also keep in mind that, intu- 
itively, the use of several clusters may introduce additional dif- 
ficulties, much like the case of a single cluster with pronounced 
substructure. 

Finally, we should stress that a potential problem for the 
application of our method could be that of multiple lensing. In 
fact, if weak lensing analysis is used (in order to solve the prob- 
lem of the sheet invariance), "tidal" effects of other clusters (or 
from large scale structures) could become important and may 
jeopardize the appUcability of the method. On the other hand, 
it is often very hard to break the sheet invariance. In view of the 
above considerations, it would be very interesting to develop a 
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method able to overcome these problems, possibly by decou- 
pling the sheet and the global scaling invariance properties. 

8. Conclusions 

In this paper we have proposed and discussed a method, based 
on weak lensing and redshift observations, to reconstruct the 
lens mass distribution and to obtain, at the same time, infor- 
mation on the geometry of the universe. Simulations based on 
the homogeneous Friedmann-Lemaitre models have shown that 
the method is viable when the number of source galaxies is 
sufficiently high. In principle, the same method could be ex- 
tended to other cosmological models, provided that the angular 
diameter-distance relation be available; some of these would be 
harder to handle because they involve a larger number of free 
parameters (e.g., the inhomogeneous "in average" Friedmann- 
Lemaitre models; see Kayser et al. 1997). 

In this paper we have used a "parameter-free" approach in 
the determination of the cosmological weight. In fact, in the re- 
construction process we have tried to determine the values of 
for a number of different redshifts (25 in our simulations), 
without any assumption on its shape. A different method could 
have been used, namely the direct determination of the cosmo- 
logical parameters without the intermediate determination of 
the cosmological weight w{z) (parametric approach). In this 
case the cosmological parameters could have been estimated 
using, for example, a maximum likelihood method. A simi- 
lar dichotomy exists between parametric and non-parametric 
mass reconstructions. We have decided to use here the non- 
parametric approach as more "conservative," but it would be 
interesting to develop a parametric method in detail for com- 
parison. 

One improvement that could be made on the current 
method is based on more direct use of the shear in the deter- 
mination of the cosmological weig ht (see Eq. (|l8|)). In fact, 
the complex shear 7(0) is not an independent quantity, in the 
sense that a relation between its partial derivatives holds (i.e., 
V A u{0) — 0; see Paper II for a detailed discussion). Hence it 
is natural to expect that the inclusion of an additional relation, 
not considered in the present paper, could improve the deter- 
mination of w{z). This might be implemented by resorting to 
a maximum likelihood approach for the measurement of the 
shear map. 
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Appendix A: Tlie slieet invariance 

In this Appendix we describe the sheet invariance for non- 
weak lenses. Our discussion largely follows Seitz & Schneider 
(1997). 

Following the approach of Appendix C, we start by tak- 
ing a fixed cosmological weight equal to the true cosmological 
weight times a given constant, i.e. w{z) = vwo{z). By using a 
global scaling invariance with scaling factor k — \/v v^t can 
make w{z) — wo{z). An approximate expression for the mean 
value of the measured shear, valid also for non-weak lenses, 
can be shown to be 



{w)^-k{w ), _ 

\i) - T~\ r^^^o 



(A.l) 



In this equation we have neglected the smoothing effect of the 
weight function W{0, 6'). From this relation we can derive the 
form of the sheet invariance. Suppose that (7) = (1 — C)7o, 
with C constant. This simply means that the measured value of 
7 is proportional to the true shear. Then the previous equation 
gives 



7Z2^ + (1 - C)k^ 



(A.2) 



This equation is consistent with the assumption 7 = (1 — C)7o 
because the "sheet constant" C{w)z/ {w'^)z does not affect the 
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shear. Equation (g5[) thus gives the form of the sheet invariance 
for the dimensionless projected density k in the general case. 

If we use these relations in Eq. (^Tj) we can derive the 
expected transformation induced on the measurement of the 
cosmological weight. The expressions obtained, rather compli- 
cated, are not reported here. However it is clear that the cos- 
mological weight changes, i.e. the measured w{z) is different 
from 'Wo{z). Moreover, it can be shown that w{z)/wo{z) is not 
constant: this of course means that the sheet invariance has a 
different effect on the cosmological weight with respect to the 
global scaling invariance. 

Appendix B: A constraint on the total projected mass of 
tlie lens 

The goal of this Appendix is to prove Eq. (following Lom- 
bardi 1996). As a first step we note that 

/ [1 - det A{e)] (fe ^ lim / [l - det A[e)\ A^Q 



lim 

R — >QO 



(B.l) 



In this equation B n is the disk of radius R centered in the origin 
of the 0-plane, while 9^{Bii) is the image, through the ray- 
tracing function 0^{0), of the disk Bj^. The second equality 
holds because of the change of variable formula. 

A well-known property of gravitational lenses is that the 
deflection j3{6) — 6^ (6) — 6 can be expressed as the gradi- 
ent of the deflection potential tp{6) (given the symmetry of the 
Jacobian matrix A of Eq. (H)). Such potential can be written as 



^(0) = i y" K{e')iii\\e -e'\\d^9' 



(B.2) 



For an isolated distribution of matter, such that the mass den- 
sity K vanishes outside a given disk Bj^, the potential can be 
expressed (outside Bjf) in terms of a multipole expansion 



where 



M = J K{e')d^e' , 



e[K{e')d^e' 



9yt 
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(B.4) 



(B.5) 



(B.6) 



(The first physically important term after M is the quadrupole 
Qij, as the dipole pi vanishes if the origin of axes is chosen in a 
suitable way.) Thus we obtain for the deflection 



1 /MO p 2{pe)e 



TT V 612 



04 



.(B.7) 



Usin g th is expression we can estimate the last integral in 
Eq. ( |B.1| ). In fact, using the definition of (3{0) and neglecting 



terms of order 0{9 ^) , we can say that the region 0^ [Bb) is a 
disk of radius 



W = R 



M 



ttR 
Thus we find 

d^e ^ ttI R 



(B.8) 



M 



0{R-^) 



ttR^ ~2M + 0(R-^) 



(B.9) 



Inserting this expression in Eq. ( |B.1| ), we obtain Eq. ([27|). A 
more detailed discussion shows that, in reality, the error in- 
volved is not 0(i?^^) but 0(^R^^), and precisely it is of order 
0{aP / R^ ,q/ R^). This property has to do with the fact that 
the dipole p, which would be responsible for the error of order 
©(i?^^), can be made to vanish with a suitable choice of the 
origin of the axes. However, the origin of coordinates is not im- 
portant when we take the limit R ^ oo, and thus the error is 
only of order 

An improved version of Eq. ( ^7| ) can be obtained by retain- 
ing the term of order O (NP/R^): 



1 r A/f^ 

- [i - det A{e)\ d^e ~ M 

2 J Br 



(B.IO) 



This expression may become useful when the quadrupole qij is 
neglegible. 

Appendix C: Expectation values and errors 

An analytic calculation of the expectation values and errors for 
the shear and for the cosmological weight is important both 
for a first check of the quality of our method and for its prac- 
tical implementation. The problem is not easy to solve in the 
general case, mainly because it involves iterations of different 
equations. For this reason we will often simplify the discussion 
by introducing some approximations. 

We expect calculations for the covariance matrix of w to 
be especially difficult, mainly because in Eq. ( p^ for w{z) the 
quantity 7(0) has been calculated using the same galaxies to 
be used for w{z). Thus the errors on w and on 7 should be 
correlated. (Strictly speaking this is true also for the covari- 
ance matrix of 7, but such dependence is not very important 
for our problem.) For simplicity, in the following we will con- 
sider either the cosmological weight or the shear map fixed in 
the calculation of the expectation value and error of the other 
quantity. The results obtained, checked by simulations, show 
that this approximation is good. 

The general method used for calculating expectation values 
and errors is described in detail in Paper I and Paper II (see es- 
pecially Appendix C of the latter article) and is based on the 
assumption that we can linearize the relevant equations near 
the mean value of the random variables. Calculations will be 
done by replacing summations with integrals, so that Poisson 
noise will be neglected. As explained in Paper II, Poisson noise 
generally adds only a small correction to the dominant sources 
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of error. For this purpose we define the redshift probability dis- 
tribution Pz{z). Hence the probability that a galaxy is observed 
in the solid angle iPd and with redshift in the range [z, z + dz] 
is given by ppz (z) (Pd dz (notice that source galaxies are taken 
to be uniformly distributed on the field, and thus p does not 
depend on 6). In the following, for simplicity, we will also as- 
sume that the distribution of source ellipticities is independent 
of redshift. 

C.l. Weak lensing 

C. 1 . 1 . Expectation values 

Calculations are not difficult in the weak lensing limit because 
all expressions are linear The shear is calculated from Eq. (p7[), 
and thus, from Eq. (|l6|), the expected value of the shear is 



(7>(^) 



E„^iW^(6),g("))w(z("))u.o(z("))7o(6)(")) 



E:=iW{e,e(-))[w{zM)Y 

If we move to a continuous distribution we find 



.(C.l) 



(7>(e) 



d^e' / dz' pz{z')w{e,e')w{z')wo{z') x 



X 7o(0') J d^e' J dz' Wi9, 0') X 

p.{z')[w{z')] 



(C.2) 



We observe now that the integrals on d^6' and on dz' factorize, 
thus leaving simply 



me) 



K 



fd^e'wie,9')-/oie') 
J d'^0'w{e,e') 



(C.3) 



The constant K is obviously related to the global scaling in- 
variance and is given by 



K 



J dz' pz{z')w{z')wo{z') 
J dz' p,iz')[wiz')f 



(C.4) 



dz' d^9'p,{z')W,{z,z') 



Jd^0"w{e',e")-fa{e") 



(C.6) 



As for the shear the integrals on the redshift factorize. Let us 
consider in more detail the integrals on the angles, i.e. 5R(- • •) 
in the numerator and j d^d' J d^9"\ ■ ■ ■ \^ in the denominator. 
These two quantities differ only by one integration with weight 
W{6', 9"). Thus we may argue that 



{w){z) 



1 jdz'p,{z')W,{z,z')wQ{z') 



{C.l) 



K J dz'pz{z')Wz{z,z') 

We note that here the constant K appears at the denominator, 
consistent with the global scaling invariance. 

In the following, in order to simplify equations, we will 
introduce the normalization condition 



dz' Pz{z')W,^{z,z') = 1 \fz 



(C.8) 



Similarly, we will consider the angular weight functions (see 
Eq. ( |C.3| )) to be normalized, i.e. to satisfy the condition 



P 



J w{e,e')d^0' ^ 1 ve 



(C.9) 



We recall that this does not limit our discussion, since only 
relative values of W are important (see Paper 11). 

C.l. 2. Errors 

Let us now consider the error (i.e. the covariance matrix) of 7. 
Using the definitions of Paper I and Paper II and calling t\ and 
£2 the real and imaginary components of the elhpticity e, we 
find 



N 



n— 1 ^ ^ ^ ^ 

In this equation we have used the matrix notation for the partial 
derivatives. For example, the last matrix, because of the trans- 
pose, has d^j{0') / de"^^ in the element with row i and col- 
umn j. From Eq. (|5]) we find that Cov(e) = Cov(e*') = cid, 
where c is a positive constant and Id is the 2x2 identity matrix 
(note that, as the source ellipticities are taken to be independent 
of redshift, c is independent of redshift as well). The partial 
derivatives of 7 can be calculated from Eq. (^tJ). By moving to 
the continuous description, after some manipulations we find 

LovH;&,& ) = CO- ^ — id, (C.ll) 

,5) J dz" p,{z")[w{z")] 



(n)^ fdW) 



.10) 



Note that K is ©-independent, so that the use of an inaccurate 
cosmological weight affects only the scale of the expected mea- 
sured shear. A result similar to Eq. ( p.3[ ) has been obtained in 
Paper II. 

From Eq. (p^), the expectation value of the cosmological 
weight is given by 



By moving to a continuous distribution and by replacing 7 with 
its expected value given by Eq. (C.3), we obtain 



{w){z) 



1 

K 



dz' d^9'p,{z')W,{z,z')wo{z')x 



I 



m d'0"w{e',9")ro{e'hoie") 



where the normalization condition ( |C.9| ) has been used. This 
result is a generalization of a similar result obtained in Paper II 
(see Eq. (26) of that article). 

If we consider 7 fixed, the covariance matrix of w is given 



N 



Cov^w; z, z') = , 
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The partial derivatives can be calculated from Eq. ([T^. Turning 
to a continuous description and using the normalization condi- 
tion dCSl) we find 



Loviw: z,z ) — 7j . (C.13) 

P fd^e"\j{e")\ 

In the case z — z', this equation has a simple interpretation. Let 
us call Nw{z) the number of galaxies in the whole field for 
which Wz(z, z^"^^^ is significantly different from zero. Then 
the variance of w{z) is given simply by c/ [Afvi'(-z)(l7p)] . 
where (|7p) is the mean value of the shear on the field. The 
covariance of u) is in general proportional to 1 / ( 1 7 P ) , and thus 
strong clusters should be preferred. However, as discussed be- 
fore, for strong clusters the sheet invariance can generate seri- 
ous problems in the determination of w{z). This simple reason- 
ing also clarifies the behaviour of Cov(w; z, z') as a function of 
the weight Wz {z, z') used. Since the expected variance of w(z) 
is proportional to 1/Nw{z), broad functions Wz{z, z') should 
be preferred. However it is clear that the coherence length of w, 
i.e. the maximum value of |z — z'| that makes Cov(w; z, z') sig- 
nificantly different from zero, is given by Nw{z): this would 
suggest the use of a function Wz with small Nw{z). 

Equation ( |C.13 ) has been calculated assuming a generic 
shear map 7(0). In practice, the shear used to calculate 'w{z) 
is the estimation 7(6), and thus 7(6) in Eq. (C.13) should be 
replaced by 7(6). 



C.1.3. Role of the global scaling invariance 

So far no assumption on w{z) has been made. If we consider 
the expectation value of 7, we see that the assumed cosmologi- 
cal weight enters through the quantity K. We have akeady ob- 
served that indeed this quantity is related to the global scaling 
invariance. In contrast to the expression for (7), the expected 
covariance of 7, given by Eq. ( |C.11| ), has a more complicated 
dependence on w{z). When K = 1 the measured {w){z) is 
closest to the true wo{z) (see Eq. ( p.7| )). From Eq. ( p.ll[ ) we 
see that the minimum of Cov(7), under the constraint K = 1, 
is attained when / dz' ■w(z')wo(z') is maximum. It is then easy 
to show, using the Cauchy-Schwarz inequality, that this occurs 
when w{z) = wo{z). 

The situation for the covariance of w is different. In fact. 



from Eq. (C.13), it is clear that the constant K does not enter 



explicitly in the expression of this covariance. This property 
will play an important role in the statistical analysis of Ap- 
pendix D. 

In conclusion, the global scaling invariance, represented 
by the quantity K, has different effects on the various quan- 
tities. This quantity enters directly in the expressions of the ex- 
pected values of 7 and w, while the expression for Cov{w) 
is not affected. However, the expected values of the shear and 
of the cosmological weight differ for an important characteris- 
tic. In fact, if the global scaling invariance could be broken, 
7 would be estimated also without the knowledge of w{z) 
(see Eq. (CJ)), while w{z) cannot be calculated without the 



knowledge of 7. Therefore, the covariance of 7, defined as 
((7(0) - (7)(0))(7(0') - (7)(0')))' can be identified with 
the eiTor {{j{9) - 70 (0)) (7(6') - 7o(0'))) even if the cos- 
mological weight used is not a good representation of wq. The 
same, of course, is not true for w, since there is no relation be- 
tween (w) and Wo, if the shear map used differs significantly 
from 70 . 

C.2. General case 

Calculations for the general case are exceedingly difficult. If 
the lens is not too strong (say, \g\ < 0.5) we might argue that 
a first approximation could be given by the results of the pre- 
vious section. A posteriori, numerical simulations (described 
in the main text) confirm the applicability of the weak lensing 
analysis to relatively strong lenses. On the other hand, if the 
lens is not weak and the sheet invariance cannot be broken, we 
have already shown that the method described here is difficult 
to apply. 



Appendix D: Extraction of cosmological information 

We now have all the tools ready, in order to proceed from the 
estimation of w{z) to the desired constraints on the cosmologi- 
cal parameters. In order to describe practically this function, we 
start by discretizing the problem in the redshift variable, which 
so far, for simplicity, has been kept continuous. We consider a 
grid of points {zi}, with i = 1, . . . ,Nz. We then introduce the 
notation 



Wi = w{zi) , 

im) = {w){z^) , 

Zij = Cov(w; Zi, Zj) , 



(D.l) 
(D.2) 
(D.3) 



for the cosmological weight function, its expectation value, and 
its covariance on such a grid. At least for the weak lensing 
limit, the expressi on f or {wj ) and Zij can now be taken to be 
known (see Eqs. {ZJ) and ( |C.13 )). A similar discretization is 



performed in the 0-space. The least squares method for the ex- 
traction of the cosmological information focuses on the study 
of the chi-square function: 



(D.4) 



with ^ = (il, k) a point in the parameter space S. The 
ambiguity associated with the global scaling invariance is rec- 
ognized and taken into account by means of the scale k which 
is treated as a free parameter in addition to the physical param- 
eters 51 and Vtiy that we would like to constrain. Note that Vl 
and JIa enter in a complicated, non-linear manner, through the 
convolution of in the expression for {w) . The study of (^) 
can be addressed with three different goals. 
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D.l. Point estimation 

Here one is interested in determining the most probable val- 
ues of the cosmological parameters. These are the values ^ that 
minimize x^, obtained by solving the set of equations 



= 



(D.5) 



This could be completed with a determination of the related 
errors. 

Here we should recall that the least squares estimation has 
the following properties (see Eadie et al. 1971): (i) If the ex- 
pectation values {wi) have been determined correctly, then the 
estimation is consistent, i.e. it leads to the true cosmological pa- 
rameters in the limit N ^ oo. (ii) If the measured weights Wi 
follow a Gaussian distribution, then the least squares method 
is equivalent to the maximum likelihood method. This occurs 
when the number of galaxies is sufficiently high (Nw{z) > 
10). (iii) If the mean {w){z) is a linear function of the pa- 
rameters, or if a linear approximation can be used, the Gauss- 
Markov theorem states that the least squares method has min- 
imum variance among the unbiased estimators that are linear 
functions of the observed quantities. This last property reas- 
sures us of the quality of the method used. Fortunately, at least 
for the scale parameter k the theorem holds. 

D.2. Interval estimation 

The interval estimation aims at identifying the relevant confi- 
dence regions in the parameter space S, starting from a given 
confidence level, CL < 1. The region is usually determined by 
referring to the chi-square probability distribution, which de- 
pends on the number of degrees of freedom. 

In the present case, given the values of the cosmological 
parameters, measurements of the cosmological weight will be 
realized (in the space of all possible measurements) with some 
probability distribution. Let us isolate in such a space a region 
with the following property: the probability that a set of mea- 
surements {wi} be realized inside it is greater than a desired 
level, the confidence level CL. Suppose that a change of vari- 
ables from {wi} to {vi} can be found such that the correspond- 
ing region in the new w-space does not depend on the cosmo- 
logical parameters. Then, given one set of measurements {wi}, 
we may associate that region in w-space to the so-called confi- 
dence region in parameter space. 

Our case requires additional discussion because one of the 
parameters, the scale k, is non-physical, and we are interested 
in defining confidence regions in the plane. In order to 

do so, we might project the region found in the S-space onto 
the plane of physical parameters. Since k can be ignored, an 
efficient approach is the following. We note that if we mini- 
mize with respect to one parameter, the new quantity fol- 
lows asymptotically (N ~> oo) a chi-square distribution with 
Nz — 1 degrees of freedom. Therefore, we may introduce the 
quantity 



which is expected to follow a chi-square distribution with two 
degrees of freedom. With this device, the confidence regions 
can be drawn in a straightforward way, independently of the 
number of points in the redshift grid. For finite values of 
N, £2 is not required to follow strictly a chi-square distribution 
with two degrees of freedom; however, £2 is always a reason- 
able statistic. 

The method described above turns out to be just an applica- 
tion of the well-known "likelihood ratio" method to our prob- 
lem (with the global scaling invariance taken into account). In 
fact, it can be shown that the quantity £2 can be expressed as 

N 
n=l 

-ln/:(6(") I f],r!A,fc,7(0(")))] . (D.7) 

In this expression C is the "likelihood" for a single galaxy (see 
Eadie et al. 1971), i.e. the probability of observing a galaxy 
with ellipticity e^"^ when the shear is 7(0) and when the "cos- 
mological" parameters are fl, Ha, and k. Then, using the cen- 
tral limit theorem and the results ( Z.l) and ( p. 13 ) for the mean 
and the covariance of the cosmological weight, the relation be- 
tween Eq. ( p.6j ) and (D.7) follows. 

The quantity £2 will prove to be useful also as a test of 
hypotheses, as explained below. 

D.3. Test of hypotheses 

The third analysis that can be performed on the measured cos- 
mological weight is the test of hypotheses. In this case a given 
hypothesis Iiq, called "null hypothesis," is tested against an al- 
ternative hypothesis hi. A test is a method to choose between 
ho and hi from a given set of observations. Hence it can be 
thought as a region in the space of the observed quantities: if 
the point representing the observations lies inside this area, the 
null hypothesis is true, otherwise it is false. 

All hypotheses that we may be interested in testing would 
actually be composite hypotheses, because of the presence of 
the scale parameter k. Let Sq be the subset of S for which 
the null hypothesis is true (for example, J7 + Oa = 1, with 
k arbitrary). Then the alternative hypothesis holds on the set 
S \ So . The likelihood ratio is defined as 



max I 







(D.8) 



£2{n,nA) 



minx iO 

k 



(D.6) 



Clearly < / < 1. Then, the likelihood ratio test is given by the 
following rule: if l < Iq then choose Hq, otherwise choose hi. 
The limiting value Iq has to be chosen so that the probability of 
"loss error" (i.e., the rejection of the null hypothesis when in 
reality ho is true) is small, and the threshold value may be set, 
for example, at 5%. Such a test is found empirically to produce 
useful and significant results. 

If we want to test a hypothesis hg fixing separately the val- 
ues of O and r^A (so that Sq reduces to a straight line), by 
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combining the discussion of Eq. ( p.7| ) with that of Eq. ( p.8| ) 
we see that we can express ^2 in terms of as 

^2 = = minx'(0 -niinx'(0 ■ (D.9) 

Thus it can be shown that, if the ho hypothesis is true, the quan- 
tity (—2 hi/) follows (asymptotically) a chi-square distribution 
with two degrees of freedom. 

For different types of ho the quantity (—2 In/) follows 
(asymptotically) a chi-square distribution with degrees of free- 
dom equal to the number of parameters fixed by ho- 

At this point, the "level of acceptance" can be obtained 
from the cumulative chi-square distribution for the relevant 
number of degrees of freedom. In particular, if we want to test 
the Einstein-de Sitter universe (51 = 1, JIa = 0) with 5% sig- 
nificance, we use the statistic 

^2=minx'(l,0,fc)-x'(e) • (D.IO) 

k 

The universe will be considered Einstein-de Sitter if £2 < 
5.991. If, instead, we want to test the flat hypothesis (Jl + fiA = 
1) with significance 5%, we use the statistic 

h - mi?^'("' 1 - ^) - xHi) ■ (D-11) 

The universe will be considered to be flat if £1 < 3.841. 

As noted earlier, when the asymptotic (N 00) limit is 
not justified, one should resort to Monte Carlo simulations in 
order to set the appropriate "level of acceptance." 



